Entanglement between quantum fields 



Daniele Teresi and Giuseppe Compagn 
Dipartimento di Fisica, via Archtrafi 36, 90123 Palermo, Italy 
(Dated: May 26, 2011) 

We analyze entanglement between quantum interacting fields. In particular, we consider the 
entanglement between the fields in the ground state of the linear a model both in its unbroken and 
spontaneously broken symmetry phases, quantified by the Renyi entropy. We find a generalization 
of the geometric entropy area law, which implies that the quantum correlations most relevant to 
the entanglement between fields are always at very short range. We find also that the degree of 
entanglement is larger in the spontaneously broken symmetry case due to the appearance, in the 
ground state, of new kinds of fluctuations. 
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Quantum entanglement, one of the most important 
characteristics of quantum mechanics, is typically ana- 
lyzed in discrete systems, e.g. spin pairs or lattices [H, 01; 
on the other hand, continuous systems, and therefore 
Quantum Field Theory (QFT), take a fundamental role 
both in Particle Physics and in Statistical Mechanics. 
Therefore, it is of interest to understand the role of en- 
tanglement in quantum continuous systems, described by 
quantum fields. Until now, the presence entanglement in 
QFT manifested itself in the phenomenon of geometric 
entropy (3n8|], i.e. entanglement between two comple- 
mentary regions of a single field. This concept was origi- 
nally introduced in the context of black- holes f3!,'3|, where 
an observer outside the horizon cannot be affected by the 
degrees of freedom of the field within. Even if the whole 
field is in a pure state, the state of the system accessible 
to the observer, that is the external degrees of freedom 
of the field, is mixed, while quantum correlations (entan- 
glement) between degrees of freedom inside and outside 
the horizon exist. The entropy of the reduced state of 
the field outside is called geometric entropy and it is a 
measure of the entanglement between the field in the two 
regions. An important property of geometric entropy is 
that it satisfies an area law 8], i.e. it is proportional 
to the surface of separation A between the subsystems; 
namely: S cx AK^ , where A is the UV cutoff introduced 
to account of the fact that the theory is assumed to be 
valid only until an energy scale A. A crucial implication 
of the area law is that the most relevant quantum cor- 
relations are localized at distances of the order 1 /A and 
this property is independent on the field being or not 
massive. 

Another QFT situation where entanglement does 
present is in the case of interacting quantum fields. Un- 
der this condition, even if the state of the whole system 
is pure, the reduced state of each field is mixed. A first 
direct consequence is that even when the system is in the 
ground state, excitation quanta of each field are present. 
As a result, if one just focuses on one of the fields, its en- 
tropy is not zero. This condition appears naturally when 
the focused field interacts with fields experimentally not 



accessible. This is, for example, the situation occurring 
in the Higgs model, with its unobserved scalar sector. 
The field entropy in such cases is a measure of the en- 
tanglement between focused and unobserved fields. This 
entanglement gives information on the quantum correla- 
tions between the fields present in the total state. 

At best of our knowledge only the entanglement be- 
tween two regions of a single field has been subject to 
investigation. Here, instead, we shall analyze the entan- 
glement between different fields, with the aim of obtain- 
ing information on the structure of the quantum correla- 
tions between the fields, and how these depend on various 
physical conditions. To this purpose we shall consider the 
ground state of the linear a model ^] both in absence and 
presence of spontaneous symmetry breaking (SSB). We 
stress that our approach differs from an effective theory 
one, because in this latter case the focused field would 
be described by an effective pure ground state, while in 
our case it is just the mixedness of its state that allows 
to get information on the entanglement, and then on the 
structure of the correlations between the fields. 

Let us preliminarily take as our system two interacting 
scalar fields cr(x) and 7r(a;), defined on the 3-1-1 dimen- 
sional spacetime, in the total ground state. The entropy 
of the subsystem cr is a measure of the entanglement be- 
tween the two fields. As a preliminary step to evaluate 
it, we require the trace of the a-th power {a G N"*") of 
the reduced density matrix of the field cr, Tr o^. Follow- 
ing a method similar to the one used in 0, [J to obtain 
geometric entropy, we shall express it as a path-integral 
on a suitable manifold. The density matrix element of 
the whole system between two configurations a' (x) , tt' {x) 
and a" {x), it" [x) is given by the following euclidean path- 
integral: 

{a"{x),TT"{x)\p\a'{x),TT'{x)) = — f PciPTre-^^f'^''^) 

Zi J 

(1) 

where Se is the euclidean action, Zi is the partition func- 
tion, and (t{x), 7r(a;) satisfy the following conditions at 
euclidean time r = 0: (7(.x, 0+) = a'{x), 7r(x, 0+) = 
7r'(af), (j{x,0^) = cf"{x), h{x,Q^) — tt"{x). In other 
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FIG. 1. Structure of the manifold on which the path-integral 
giving Tr is performed. Planes are represented by contin- 
uous black lines, sheets by dotted red ones. 

words, the path- integral is defined in a spacetime with 
a cut along r = 0. On the two edges of the cut the 
fields coincide with the configurations (t"(x), 7r"(x) and 
(t'{x), tt'{x). In order to obtain Tr p^, first we must trace 
out the field tt; this can be done by taking identical values 
for this field on the two edges of the cut and then integrat- 
ing upon all its possible boundary configurations. This 
amounts to consider a{x) as living on a spacetime with 
a cut at T = 0, whereas t:{x) lives on ordinary spacetime 
without any cut. Within this approach {cr"{x)\ \<j'{x)) 
is represented by a path-integral on a manifold made up 
of two of these spacetimes: <j{x) lives on the manifold 
obtained sewing together two opposite edges of the cuts 
belonging to the two different spacetimes. On the re- 
maining two open edges u{x) is bound to the configura- 
tions u' {x) and a"{x). The manifold for 7r(a;) is obtained 
simply duplicating the ordinary spacetime, without any 
cut or bond. Finally, Tr is obtained sewing together 
the remaining open edges and thus the manifold obtained 
is given by the structure sketched in Fig.[T] The manifold 
on which the path-integral is calculated, not containing 
conical singularities, is simpler than that used to obtain 
geometric entropy [5, X]- This allows to treat our interact- 
ing field model. Since the path-integral for is obtained 
using the product of two expressions of the form ([l} , Zi 
will appear squared, and the final result is that: 

Trp,^ = ^|paP.e-^-(--) = ^ (2) 

where the path-integral Z2 is calculated on the structure 
of Fig. [TJ Generalization to Trp^, a £ is trivial: 
Trp^ — Za/Zi, where Za is the path integral on the 
generalization to a > 2 of the structure in Fig. [TJ If the 
system contains more than two interacting fields and we 
want to focus on one, thus analyzing the entanglement 
between this field and the rest of the system, the remain- 
ing fields must be traced out. The manifold is like tt's 
one for traced out fields and like cr's one for the focused 
field. In the following, we shall first study a two-field 
model with cubic interactions and then generalize to the 
linear a model with N >2 fields. 

The manifold for the functional integration of the tt 
field consists of a disconnected spacetimes, which will 
be referred to as planes, that for a, instead, of a dis- 
connected pieces named sheets, which have the struc- 
ture of ordinary spacetimes. A single plane is an or- 
dinary euclidean spacetime, and therefore the propaga- 




FIG. 2. O(A^) dia grams that contribute to Sa- cr's propa- 
gators are denoted by continuous lines, tt's ones by dashed 
lines. 

tor inside a given plane is the ordinary euclidean one: 
K{x,y) — J ^'^ k ^^^'i Q^K^-v) ^ where m is the mass of 
the field. Because planes are disconnected, instead, the 
propagator of tt between two different planes is zero. 
Analogously, as concerns a, the propagator D{x,y) in- 
side a sheet is again the ordinary one, whereas it is zero 
between different sheets. 

In our two-field model, let us take the interaction 
terms as given by the following euclidean interaction La- 
grangian: Cj = Xva^ + Xvan'^, where v has the dimen- 
sions of a mass. For sufficiently small coupling constant 
A, it is possible to expand the exponential in ([2]) in Taylor 
series, and express every term by vacuum Feynman dia- 
grams with (7(7(7 and amr vertices. Naming Zi and Z^ 
the path integrals Zi and Za in the free theory, and intro- 
ducing Wi , Wa as follows: Zi = Z°e^' , Z^ = Z^e^- , 
only connected diagrams will contribute to Wi and Wa- 
Z^ factorizes in the free path-integrals on 17 and on tt, 
which in turn are the product of the contributions from 
every single plane or sheet. Since in the free theory 
the path integral on a single plane/ sheet is the usual 
one, ZO = {Z^y, and therefore: Trp^" = gW^-aW^_ 
Tr p" can be directly used to obtain the Renyi en- 
tropy of the focused field a: Sa — j^lnTrp^ = 
^^(a Wi — Wa); this is a generalized entropy that gives 
a measure of the mixture of the state and therefore of the 
entanglement between the subsystems. 

The first contribution to Ws, and thus to Sa, is 
O(A^); at this order only the first diagram in Fig. [2] 
contributes. The other O(A^) connected diagrams (not 
shown in Fig. [J) have two vertices connected only by 
(7 propagators and their contribution to Wa is equal to 
a times the contribution to Wi. Evaluating Wa, for the 
first vertex in a given sheet, the integration on the second 
must occur on the same sheet and then all the tt and a 
propagators are standard. The contribution on this sheet 
is equal to Wi 's one. Because there are a sheets, the con- 
tributions of each of these diagrams to Wa and aWi are 
equal, and cancel in Sa. This argument confirms the 
physical expectation that only virtual processes in which 
both the fields propagate between the two vertices can 
contribute to the entanglement between the fields. 

Let us derive the contribution given by the first dia- 
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gram in Fig. [2] to Wa- Because this diagram is connected 
both by a and tt propagators, the integration of the sec- 
ond vertex must be performed, as said before, both on 
the same sheet and plane of the first. For the first vertex 
on a given plane with r > (r < 0), the second must 
then be integrated also in the same region t > (r < 0), 
and the propagators are standard. The first vertex can 
be on any of the a planes, thus the integration on the two 
vertices must be performed on a couples of concordant 
semi-spacetimes. In Wi these integration are performed 
on a single full spacetime, and therefore in a Wi — Wa 
only the integration on the a couples of discordant semi- 
spacetimes survives, and Sa up to O(A^) is: 
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^T>0 T<0 T<0 T>0 ^ 

D{x,y)K{x,y)K{x,y) (3) 



We will apply the above technique and results to the 
linear a model, both in its unbroken and spontaneously 
broken symmetry phases. 

Unbroken symmetry. The linear a model consists of 
N scalar fields with a quartic interaction, and has 0{N) 
symmetry. In the unbroken symmetry case, naming a 
the focused field and nk, k = 1, . . ._jjy — 1 the remaining 



ones, the euclidean Lagrangian is 
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where is the counterterm Lagrangian. All the fields 
have mass to, and each vertex has 4 lines. Similarly to 
the case of cubic interactions, only the second diagram 
in Fig. [2] contributes at order A^ [l2| . The expression for 
this diagram is similar to ([3|), but with an extra D{x,y) 
(caused by the extra a line) and a coefficient {N — l)A^/2 
instead of X^v^; the factor A^— 1 being due to the presence 
of the — 1 fields tt'' . Expressing the propagators in the 
momentum space, integrating on the positions x, y of the 
vertices and on the total momentum, we get: 
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N ~1 aX^ 



-V 



d^pd^kdH- 
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2 a — I JJJ ^ k'^ + P + m? 

1 {p + k + If - {po + ka + kf + m'^ 



P'' + TO'^ 



{p + k + iy + to2 [{p + k + If + to2 



where the subscript denotes the time component, V is 
the 3D volume and the last factor is connected to the 
space integration on discordant semi-spacetimes. Using 
Pauli-Villars regulator method jUj] and keeping only the 
dominant A divergent term, we get: 
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(4) 



The prefactor C is a finite adimensional positive number 
whose precise value, similarly to what happens for ge- 
ometric entropy in the multi-dimensional case, depends 
on the details of the regularization procedure, that is on 
the unknown details of the theory at energies of the or- 
der of (or larger than) the cutoff A, and therefore it is 
not an universal quantity (in the sense of statistical me- 
chanics). In spite of this, it will be yet possible to draw 
some physical consequences from s'^\ in particular on 
the structure of the quantum correlations between the a 
and TTS fields. It is possible to show that, if the mutual 
information is used instead of the entropy to quantify the 
entanglement [l^], a similar not universal prefactor ap- 
pears. Moreover, one may expect that, analogously to the 
case of geometric entropy [IJ], the neglected terms con- 
tain finite terms and/or logarithmic divergences whose 
prefactors are insensible to the regularization procedure, 
and thus universal. 

Before going ahead to examine the physical conse- 
quences of (HI, we observe that Von Neumann entropy 
S is typically used to study the entanglement proper- 
ties of fields. For geometric entropy <S can usually be 
obtained from the Renyi entropy Sa using the so-called 
replica trick 0-0] • This amounts to analytically continue 
Sa from integer to real a and obtain S by the relation 
limQ,_).i Sa = S. If we tried to use this method by ana- 
lytically continuing ^ to real a in the most natural way 
(promoting a to real values), and take the limit a — > 1, 
we would obtain infinite for S. This does not necessar- 
ily mean that S is truly infinite, but it may be due to 
the fact that analytic continuation is not uniquely de- 
fined if one only knows Sa for integer ot. In particular, 
that this trivial analytical continuation is not the correct 
one is confirmed by the fact that the equivalent relation 
S — — lima -1.1 ^ Tr would give an S different from 
the previous one and smaller than 52, whereas it must 
be always S > S2. This shows that the replica trick can- 
not be used in our case. Von Neumann entropy cannot 
be obtained, and we must use directly the Renyi entropy 
(with integer a). 

Now, we shall examine what physical implications can 
be extracted from Q. First, we note that S'^^ is propor- 
tional to the number of tt'^ fields; this is due to the fact 
that at this order the virtual processes that contribute to 
it involve the interaction between a and the tt*^ each at 
a time. Therefore, being each of the A" — 1 tt'^ involved 
only in an independent process, their contributions sum. 
Second, proportionality to Vh^ means that the entropy 
is extensive and that the most important correlations to 
the entanglement between a and tt*^ fields are those at 
a distance of the order of 1/A. In fact, Vh^ is of the 
order of the number of degrees of freedom of the field tr, 
and (HI says that each of these contributes independently 
to the entropy, without any other length scales different 
entering. This does not mean that quantum correlations 
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with greater length scale, for example of the order of 
1/m, are absent, but that they are small compared to 
those at length scale 1/A. Finally, because of the 0{N) 
symmetry of the theory, the same result would have been 
obtained if we had focused on any tt field, tracing out the 
remaining tts and a. 

Spontaneously broken symmetry. In the SSB case 
— < 0, the a field acquires the vacuum expec- 
tation value V = and it has to be expanded around 
this value; formally this can be obtained by the shift 
a ^ v + a. After this, cr and tt*^ are the fluctuating fields 
around the classical vacuum, and the Lagrangian takes 
the form tnl: 



the physical ones. We thus finally obtain: 
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In this case naive perturbative expansion confronts two 
problems. First, the presence of the massless Goldstone 
bosons TT*^ cause higher and higher infrared divergences 
with the increase of the perturbation order. Second, 
there are tadpole subdiagrams whose tail, carrying zero 
momentum, does not lower the degree of divergence and 
which are not canceled by the counterterms. These sub- 
diagrams cause the appearance of higher and higher UV 
divergences. To avoid the last problem, we should resum 
the contributions of all the possible insertions of tadpoles 
in a given diagram. An equivalent, but simpler, way to 
obtain this is to expand the field around a value u ^ v 
such that the whole contribution of the tadpole insertions 
becomes finite or even null, and a perturbative approach 
is possible. We are free to choose a shift value u ^ v 
because entanglement is a property of the structure of 
the state and it does not depend on how we name the 
states; shifting a field by a constant is equivalent to label 
in a different way the states of the basis of eigenvectors 
of the field operators, without mixing the subspaces of 
the subsystems, and this does not change the entropy 
of the subsystem. Following this approach, explicit cal- 
culations show that the valuse that makes possible the 
perturbative approach is ~ Ct , where Ct 

is a positive constant and A„ is the physical coupling 
constant at energy u in the standard MS renormaliza- 
tion scheme 11 1| . with the assumption Au ^ 1. Expand- 



ing around u, a and tt*^ respectively acquire the formal 
masses ~ 3Au^, ~ Xu^ ^ 0. The appearance of a 
mass for the tt fields solves the IR divergences problem. 
Now both diagrams of Fig. [2] contribute, at the lowest 
order, to Sa- These contributions can be obtained in a 
way similar to the unbroken symmetry case, using u in- 
stead of w, and the formal masses mt and rfii. instead of 
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N + 8 lnl/A„ 



(5) 

where Cb is a positive constant with the same character- 
istics of C. The first contribution to the SSB entropy is 
the same of the unbroken symmetry's one and thus en- 
tanglement in the SSB case is increased. The reason is 
that new kinds of fluctuations, involving the background 
field that breaks the symmetry, are present and the sec- 
ond term in ^ represents the leading contribution. It is 
non-analytic in A„ and this non-analyticity causes prob- 
lems in the naive perturbative expansion. At variance 
with the unbroken symmetry case, the dependence on N 
is not simply a factor iV — 1 . The reason is that by resum- 
mation we are taking into account processes at arbitrary 
high order, included those in which there are interactions 
between different tts that, therefore, do not contribute in- 
dependently anymore. The resolution of the problems of 
the naive perturbation expansion has given rise to the 
appearance, in place of the length-scale of a much 
shorter length-scale 1/u of the order of A~^ In^^^ 1/A„. 
Finally, at variance with the unbroken symmetry case, 
the reduced entropy of a single tt*"', obtained tracing out 
a and the remaining tts, would be different from ([S]), be- 
cause the 0{N) symmetry is broken. At this order the 



change is simply a factor instead of 



(jV-l)^ 



the 

second term of ([S]). 

Conclusions. We have studied entanglement between 
interacting quantum fields in the case of the linear a 
model in the ground state, both in its unbroken and 
spontaneously broken symmetry phases. Entanglement, 
as measured by the Renyi entropy, has been shown to be 
larger in the SSB phase, and results proportional to the 
volume in both phases. These results, although they con- 
tain not universal prefactors, give information about the 
physical structure of the quantum correlations between 
the fields. Proportionality to the volume can be inter- 
preted as a generalization of the geometric entropy area 
law, with the volume being a sort of generalized separa- 
tion surface between the subsystems, i.e. the fields. This 
implies, similarly to the geometric entropy case, that the 
most relevant correlations are at very short range (of the 
order of A^^). In fact, they contribute within volumes of 
the order of 1/A^, whose number increases as the area, 
for geometric entropy separation surfaces, or as the vol- 
ume, in our case in which the subsystems are the fields all 
over the space. Remarkably, this result holds too in the 
SSB case, where long-range interacting massless Gold- 
stone bosons are present. Finally, the increasing of the 
entanglement with SSB is linked to the presence of new 
kind of fiuctuations, other than those of the unbroken 
symmetry phase. Although entropy of entanglement is 
difficult to measure even in discrete systems, it can be 
in principle computed by lattice numerical simulations. 
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that would presumably confirm our prediction that the 
correlations most relevant to the entanglement are those 
at very short range. 

Because the entropy of a field, due to its entanglement 
with other fields, would be sensitive to all fields in the 
theory, even those which have not yet been directly ob- 
served, it would be interesting to extend the above anal- 
ysis to more realistic cases such as the Higgs model. 
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